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ABSTRACT 

We analyze the finite temperature phase diagram of QCD with fermions in the 
adjoint representation. The simulations performed with four dynamical Majorana 
fermions show that the deconfinement and chiral phase transitions occur at two 
distinct temperatures. While the deconfinement transition is first order at Td we 
find evidence for a continuous chiral transition at a higher temperature T c ~ 8 Td. 
We observe a rapid change of bulk thermodynamic observables at Td which reflects 
the increase in the number of degrees of freedom. However, these show little vari- 
ation at T c , where the fermion condensate vanishes. We also analyze the potential 
between static fundamental and adjoint charges in all three phases and extract the 
corresponding screening masses above T d . 



1 Introduction 



The interplay between confinement and chiral symmetry restoration is one of the 
most puzzling problems in finite temperature QCD. Both phenomena seem to refer 
to different non perturbative mechanisms, where a priori unrelated length scales are 
involved. But speculations, that QCD would therefore undergo two distinct phase 
transitions [0] - deconfinement of quarks and gluons at followed by the restoration 
of chiral symmetry at some higher temperature T c - could not be confirmed by lattice 
calculations. The two effects seem to be tightly coupled in QCD. At a unique 
critical temperature the chiral condensate vanishes and the asymptotic value of the 
heavy-quark potential drops rapidly indicating deconfinement. 

Thermodynamic quantities calculated in the vicinity of the QCD phase transition 
temperature exhibit the deconfining as well as chiral symmetry restoring aspects of 
this transition. Observables like the energy or entropy density show a sudden rise 
at T c which reflects the liberation of many new partonic degrees of freedom due to 
deconfinement. In this respect the transition is very similar to that in the pure gauge 
sector where quarks are absent and chiral symmetry thus does not play a role. In 
QCD with nf light quark flavours, however, also the chiral condensate drops rapidly 
and vanishes at T c for vanishing quark masses. In fact, it generally is believed that 
details of the QCD phase transition, e.g. the order of the transition, is controlled 
by chiral symmetry. At least qualitatively this has been confirmed by numerical 
calculations. While the transition is first order for QCD with three or more light 
flavours it does seem to be continuous for two flavours as expected from general 
universality arguments @ , although at present the analysis of critical exponents for 
two flavour (n/ = 2) QCD is not entirely in agreement with expectations based on 
the universality class of 3-dimensional 0(4) spin models f|, [5[]. Also characteristic 
chiral singularities, which are expected to show up in the quark mass dependence of 
the chiral susceptibility for T<T C || have so far not been observed. 

In how far are the chiral properties of the transition dependent on the confinement- 
deconfinement transition and in how far is chiral symmetry restoration reflected in 
the behaviour of bulk thermodynamic observables? In order to get further insight 
into these questions one may analyze QCD related models in which both transi- 
tions fall apart. The SU(N) gauge theories with fermions in the adjoint representa- 
tion (aQCD) offer this possibility. Early investigations of SU(2) with twoQ adjoint 
fermions |7|, |J indeed suggest that the two phase transitions are well separated with 
Td < T c . This may not be too surprising as aQCD possesses two distinct global sym- 
metries. Unlike fermions in the fundamental representation the adjoint fermions do 
not break the global Z(N) center symmetry of the SU(N) gauge group. The model 

a In the following we will adopt the QCD convention when counting the number of fermion 
flavours, i.e. we count the number of Dirac fermion species in the adjoint representation rather 
than Majorana fermions. 
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thus has the Z(N) center symmetry and a SU{2nj) chiral symmetry || At low 
temperatures the latter is broken down to SO(2rif) while the former is expected to 
become broken at high temperature. Related to these symmetries and their break- 
ing or restoration are two order parameters, the Polyakov loop expectation value 
and the expectation value of the fermion condensate, which can characterize the 
different phases of the theory^. 

In this paper we concentrate on an analysis of the thermodynamics of aQCD. We 
determine the critical temperatures, analyze the order of the transitions and give an 
estimate of the latent heat at the deconfmement transition. Particular emphasis is 
put on the investigation of the thermodynamics in the intermediate phase, i.e. for 
Td < T < T c , where we analyze the temperature and fermion mass dependence of the 
fermion condensate, its susceptibility as well as potentials for static fermions in the 
fundamental and adjoint representations. Before discussing the thermodynamics of 
aQCD in section 3 we will introduce the model and give a description of numerical 
aspects of our calculations in section 2. Section 4 contains our conclusions. 

2 SU(3) gauge theory with fermions in the adjoint 
representation: aQCD 

In this section we want to introduce the basic properties of the lattice formulation 
of the SU(3) gauge theory with fermions in the adjoint representation. We will 
use for our investigations the standard Wilson one-plaquette gauge action. In the 
fermion sector we use the staggered fermion formulation. The only modification 
compared to QCD calculations is that the fermions are in the adjoint (8-dimensional) 
representation of the SU(3) group. As they carry now eight colour degrees of freedom 
we also need an 8-dimensional representation for the gluon fields in the discretization 
of the fermionic action, 

S = (3S G + Y,^M(U^) x ^ y , (1) 

where M(U^ 8 ') denotes the standard staggered fermion matrix with the 3-dimensional 
gauge field matrices replaced by the 8-dimensional representation U^ 8 \ (3 = 6/g 2 
is related to the inverse coupling of the gauge fields and Sq denotes the gluonic part 

b We note that the model we are discussing here is closely related to supersymmetric gauge 
theories. For a recent attempt to simulate supersymmetric Yang Mills theories on the lattice 
see [pd| . The possibility of the breaking of the Z(N) center symmetry in finite temperature 
supersymmetric models has recently also been discussed in Ref. |l2| . 
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of the action, 



S G = X:{l-^eTrf/(?)} . (2) 

We note that the 8- dimensional representation may be represented in terms of the 
3- dimensional one as 

U { 2 = Itr^X^} , (3) 

with A a being the Gell-Mann matrices and tr the 3-trace. This also shows explicitly 
that the action for adjoint fermions does not break the Z(3) center symmetry. In the 
fermionic part of the action the matrices £7® appear only in complex conjugate pairs 
Jj( 3 \ . The two center elements appearing in a Z{3) transformation therefore 
cancel each other. Like in the 577(3) gauge theory the Polyakov loop, 



1 1 Nr 

a x a x xo=l 

thus still is an order parameter of the theory for the spontaneous breaking of the 
center symmetry. As its value is related to the asymptotic behaviour of the Polyakov 
loop correlation function, 

eM-Vs(x,T)/T) = {L 3 (6)Lt(x)) ^ L\ , (5) 

which defines the potential between static fundamental charges, it is obvious that the 
thermal medium generated by aQCD will be confining for fundamental charges at 
low temperatures as long as the Z(3) center symmetry is not spontaneously broken. 
In contrast to this the expectation value of the adjoint Polyakov loop[], 

1 1 Nt 

c X c X Xo=l 

is non-zero for all values of j3. Adjoint charges are thus screened and the correspond- 
ing potential reaches a constant value at large distances^, 

exp(-V 8 (f,T)/T) = (L g (0) L 8 (x)) ^ L\ . (7) 

c To evaluate adjoint traces we may use the relation TrU^ = \TtU^\ 2 — 1, which can easily be 
derived from equation |^. 

d Note that this screening is not only related to the fact that adjoint charges appear now as 
dynamical degrees of freedom in aQCD. Even in the quenched case, i.e. the pure SU(3) gauge 
theory ad join t charges are screened at low temperatures. String breaking has been analyzed in 
this limit g|. 
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We note that the matrices are real and so is the entire fermion matrix M and its 
determinant. As a consequence of this pseudo fermion fields $, which are generally 
used to represent the fermion determinant as a Gaussian integral over bosonic fields, 
can be chosen real. This reduces the number of fermion doublers by a factor of 2, 
i.e. we can set up an exact fermion algorithm for two adjoint fermion flavours (see 
footnote a). 

The pseudo fermion action, used in our simulations thus reads 

S = f3j2h-\^TiU^) + ^ t {M(U {s) ) t M(U {8) )y 1 ^ . (8) 

For our calculations we use the hybrid Monte Carlo algorithm. We note that the 
entire molecular dynamics evolution of the system can be expressed in terms of the 
3-dimensional matrices and the fermion fields. An explicit representation of the 
matrices is only needed in the conjugate gradient routine used for the inversion 
of the fermion matrix. Technically the largest problem this prescription causes is 
the evaluation of the time derivative of M t M, which is needed in the molecular 
dynamics steps. As the fermion matrix now is quadratic in time derivatives 
generate also twice as many terms in the Hamiltonian evolution equations as in 
ordinary QCD simulations. 

Our simulations were performed on lattices of size x iV T with N T = 4 and N a = 8 
and 16. On the smaller lattice we performed simulations for 15 /^-values in the 
interval [5.2; 6.5] each for the bare masses m = 0.02, 0.04, 0.08, 0.10. In some cases 
additional runs were performed at a smaller fermion mass, m = 0.01, as well as at 
intermediate masses, m = 0.03 and m = 0.05. On the larger lattice calculations were 
performed with the fermion mass m = 0.02 for 8 /3-values between 5.25 and 6.5. At 
some selected values of the gauge coupling, /3 = 5.25, 5.4, 6.2, 6.5, which were chosen 
in the three different phases of thermal aQCD, we also calculated the fundamental 
and adjoint Polyakov loop correlation functions and extracted the corresponding 
potentials from them. 

We used the hybrid-$ algorithm || with the standard gluon and staggered fermion 
action only modified for the adjoint representation, and the conjugate gradient 
method for the inversion of the fermion matrix. The residuum was 10~ 8 inside 
the molecular dynamic steps and 10~ 13 elsewhere. The length of a trajectory has 
been fixed to At = 0.25. The number of subdivisions of a trajectory has been chosen 
such that we reached an acceptance rate of 60 — 80% which could drop to 50% in 
the vicinity of [3&. It turned out that the number of subdivisions needed mainly 
depends on the fermion mass and shows little dependence on (3. More details on our 
simulation parameters are given in Table |l| 
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mass 


#mol.dyn. steps 


#traj. for therm. 


#traj. for measurement 


8 3 x 4 lattice 


0.01 


64 


500 


1000 - 3000 


0.02 


40 


500 


1500 - 6500 


0.03 


32 


500 


1500 


0.04 


28 


500 


2800 - 7500 


0.05 


25 


500 


2500 - 6500 


0.08 


21 


500 - 1000 


1500 - 7000 


0.10 


20 


500 - 1000 


400 - 14000 


16 3 x 4 lattice 


0.02 


50-75 


100 - 1000 


500 - 2300 



Table 1: Parameters used in the hybrid Monte Carlo calculations on different lattice 
sizes and for different values of the fermion mass. 

Plaquette and Polyakov loop expectation values have been calculated on every con- 
figuration, the fermion condensate only on every 5th for the small lattice and every 
second configuration for the larger lattice. On these configurations we also calculated 
the disconnected part of the fermionic susceptibility M, 



Xr. 



N*N T 



((TrM -1 ) 2 ) - (TrM -1 ) 2 



(9) 



for which we used a noisy estimator [13| with 25 random source vectors. 



Correlation functions of Polyakov loops for fundamental and adjoint static fermions 
have been calculated in order to extract the potential for both types of fermions. The 
correlators were measured every second iteration at distances R = \jx\ + x\ + x\, 
Xi = 0,1,2,3,4 and also at on-axis distances 5, 6, 7 and 8. After averaging over 
correlations at equal distance we obtained the potentials at 33 points for 1.0 < R < 
8.0. 



3 Thermodynamics of aQCD 

As discussed in the previous section the Polyakov loop, L 3 , is an order parameter 
for the spontaneous breaking of the Z(3) center symmetry of aQCD and easily is 
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related to the confining/deconfining properties of the static potential for fundamen- 
tal charges (Eq. Less obvious is that this also reflects the deconfinement of the 
dynamical degrees of freedom of aQCD, i.e. the adjoint fermions. In order to ana- 
lyze this question we have calculated various observables reflecting thermodynamic 
properties of the hot medium which we will discuss in the following. We will start 
with a discussion of signals for the deconfinement and chiral transitions and then 
discuss thermodynamic observables related to the equation of state of aQCD. 



3.1 Deconfinement phase transition 




0.5 




5.2 5.4 5.6 5.8 6.0 6.2 6.4 P 6.6 

Figure 1: Expectation value of the absolute value of the Polyakov loop for the 
fundamental representation, L 3 , versus coupling, (3 = 6/g 2 , calculated on lattices of 
size 8 3 x 4 (a). Shown are results from simulations with different fermion masses 
and a result from simulations in the pure gauge sector (m — > oo). Figure lb gives a 
comparison of results from two different lattice sizes for m = 0.02 and in addition 
also shows the adjoint Polyakov loop Lg. 

In Figure [l] we show results for the Polyakov loop expectation values in the funda- 
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m 


0.10 


0.08 


0.04 


0.02 


0.0 




5.342 (2) 


5.312 (2) 


5.279 (2) 


5.256 (2) 


5.236 (3) 



Table 2: Critical couplings for the deconfinement transition estimated on lattices of 
size 8 3 x 4 and for different fermion masses. The last column shows the result of a 
linear extrapolation to the zero fermion mass limit. 

mental and adjoint representation. The Polyakov loop for fermions in the fundamen- 
tal representation shows a clear jump around (3 = 5.3. The size of the discontinuity 
does not indicate a significant dependence on the fermion mass in the range covered 
by our simulations. It is, however, much larger than in the quenched (m — ► oo) 
limit, which also is shown in Figure [I]. This is in accordance with the stronger 
screening observed for the static heavy fermion potential (see discussion below). We 
also note that the adjoint Polyakov loop changes from small but non-zero values to 
large non-zero values above T<j, which also signals a sudden change in the potential 
for static adjoint charges. 

Using the reweighting method of Ferrenberg-Swendsen [14[] we have located the peak 
of the Polyakov loop susceptibilities, 

XL = Nl((j-Y J L,{x) 2 )-Li\ , (10) 

at fermion masses m = 0.02, 0.04, 0.08, 0.10. This gives estimates for the critical 
couplings^ which were used for an extrapolation to the m = limit with a general 
ansatz, f3&{m) = /?d(0) + c-m a . It turns out that the exponent a is unity within sta- 
tistical errors. Finally, we thus performed an extrapolation using a = 1. This yields 
Ai(0) = 5.236(3). We will use this value of the critical coupling in the next subsec- 
tion to determine the relation between the critical temperatures for deconfinement 
and chiral symmetry restoration. 

In Figure § we show the potentials for static charges in the fundamental and adjoint 
representations calculated at two temperatures below and above Td on lattices of 
size 16 3 x 4. They clearly show the sudden change from a strictly confining (funda- 
mental charges) or confinement + string breaking (adjoint charges) potential at low 
temperatures to a screened potential above T d . We note that for T<T d the breaking 
of the adjoint string does set in quite early, i.e. already at distances R ~ 0.6/T. For 
smaller distances the linear rise of the fundamental and adjoint potentials is com- 
patible, i.e. we do not see indications for a Casimir scaling as it has been observed in 
pure SU(3) gauge theories at zero temperature [JTJ]. Below T d the potentials shown 

°We did not aim at a more precise determination of the critical points and thus did not analyze 
further their volume dependence. 
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a 


fi/T 


fundamental 


0.11 (4) 


6.4 (6) 


adjoint 


0.28 (9) 


5.6 (6) 



Table 3: Screening masses in units of the temperature, fi/T, and the Coulomb 
coefficient of the static potentials for fundamental and adjoint fermions at /3 — 5.4 
calculated on 16 3 x 4 lattice with adjoint fermions of mass m = 0.02. 

in Figure |2| have, in fact, been calculated at a temperature quite close to T d . This 
can be estimated by using the 2-loop /9-function for aQCD, 

with coefficients bo, b\ appropriate for SU(3) with two flavours in the adjoint rep- 
resentation, i.e. b = 3/(167r 2 ) and b\ = 115/(3847r 4 ). From this we conclude that 
our calculation at (5 — 5.25 for m = 0.02 corresponds to a temperature T ~ 0.97T d . 
At this temperature we still find quite a large value of the string tension. A fit to 
the potential for fundamental charges yields y/a(T)/T = 1.98 (4), which should be 
compared to a calculation in the pure gauge sector close to the critical temperature 

HH which gave ( y/a(T)/T\ = 0.93 (3). This shows that the deconfinement 

V /SU(3) 

transition in aQCD is strongly first order with a large discontinuity in the string 
tension. We will find further support for this conclusion from the analysis of the 
latent heat at Td which will be discussed in section 3.3. 

Above Td the static potentials are strongly screened. We have tried to extract a 
screening mass by fitting the potentials with an ansatz for a Debye screened Coulomb 
potential, 

m = _JL (12) 
T RT v ' 

The results from such a fit are given in Table ^. Although the statistical errors are 
of the order of 10% one sees that the screening masses for fundamental and adjoint 
charges are in agreement with each other, while the coupling strength is different. 
We note that similar studies in the pure gauge sector gave screening masses which 
are about half as large, \ijT ~ 2 — 3. This gives support to a picture where screening 
of the fundamental and adjoint potentials in aQCD results from the exchange of 
gluons with an effective thermal mass that is larger than in QCD because it now 
receives contributions from internal gluon and adjoint fermion loops. 



9 



16 
15 
14 
13 
12 
11 
10 
9 
8 



V(R)/T, T<T d 
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fundamental 
adjoint • 



1.0 2.0 3.0 4.0 5.0 6.0 



0.00 
-0.05 
-0.10 
-0.15 
-0.20 
-0.25 
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V(R)/T, T d <T<T c 



(b) 



fundamental 
adjoint 



1.0 2.0 



3.0 



4.0 



5.0 R 6.0 



Figure 2: Static potentials for fundamental and adjoint charges calculated below 
(a) and above (b) on lattices of size 16 3 x 4 at (3 = 5.25 and 5.4, respectively. 
Calculations have been performed with fermion mass m = 0.02. The lines show fits 
as discussed in the text. For (3 = 5.25 the potentials have been normalized at R = 1. 
Also note the different scales in (a) and (b). 



3.2 Chiral phase transition 



Let us now turn to a discussion of the fermion condensate (Figure |3]). It clearly feels 
the deconfinement transition and is discontinuous at j3&. For ft > @ d it, however, 
remains nonzero and large. For non-zero fermion masses this is, of course, expected. 
In order to decide, if the chiral symmetry indeed remains broken above Td, one has 
to extrapolate to the chiral limit, m — 0. However, before doing so we want to 
discuss the fermion mass dependence of the fermionic susceptibility, Eq. |]. At the 
critical temperature of a second order chiral phase transition the susceptibility is 
expected to diverge like m _ ^ 1_1// ' 5 ' ) , where universality relates the critical exponent S 
to that of the 0(6) model in three dimensions (5 ~ 5). For finite fermion masses this 
singular behaviour shows up as a peak in the susceptibility at some pseudo-critical 
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coupling f3 c (m) at which the peak height shows the same fermion mass dependence, 
Xm.max ~ m _ ( 1-1 '^. At least qualitatively this behaviour, i.e. pronounced peaks 
of the susceptibility at some pseudo-critical coupling, has been observed in QCD 
with light quarks in the fundamental representation. The close relation to the a- 
models in three dimensions established by universality arguments also suggests that 
the fermionic susceptibility diverges below T c in the zero fermion mass limit like 
mT x l 2 ||. For small fermion masses the divergence at T c will, however, be stronger 
(1 — 1/5 ~ 0.8) and the chiral transition should still be signaled by a peak in 
the susceptibility. So far, the presence of the divergence of Xm below T c could 
not be verified in QCD. To some extend this may be due to the strong shift of the 
pseudo-critical couplings with quark mass, /3 c (m), which makes an analysis of Xm for 
constant temperatures below but close to T c difficult in QCD. As will become clear 
below these characteristic singularities induced by the fluctuations of the Goldstone 
modes are, however, clearly visible in aQCD just because of the clear separation of 
the deconfinement and chiral phase transitions. 

In Figure |3] we show the fermionic susceptibilities of aQCD for several values of the 
fermion mass. First of all it is evident that the singularity related to the breaking 
of the Z(3) symmetry is also clearly visible in Xm- For temperatures above the 
susceptibility does, however, show quite an unusual behaviour not easily visible in 
QCD simulations with fermions in the fundamental representation. Above there 
exists a long, strongly mass dependent plateau at the end of which a second peak 
starts developing for the smallest fermion masses used in our simulation, i.e. for 
m = 0.02 and m = 0.01. As can be seen in Figure ^ this plateau indeed rises like 
1 / yfm as expected from the analysis of Goldstone modes in the three dimensional 
cr-models JJ. Only close to T c the susceptibilities start rising faster and thus reflect 
the more singular behaviour at T c . In this region one thus would expect that the 
appropriate form to extrapolate the fermion condensate is given by the ansatz, 

= ao + a im 1/2 . (13) 

This is indeed supported by the numerical results shown in Figure |3|. For larger 
fermion masses, of course higher order corrections become relevant. In the inter- 
mediate regime 5.4 < (5 < 5.8 we therefore have analyzed the data for the fermion 
condensate using the general ansatz 

( , 0'i/>)(/3, m) = do + aim 1 / 2 + aim + a 2 m 2 + a 3 m 3 . (14) 

For couplings (3 > 5.8 we have used this ansatz with a\ = 0. The resulting extrapola- 
tion for the fermion condensate is also shown in Figure |^. This clearly demonstrates 
that the fermion condensate remains non-zero above and only vanishes at a larger 
temperature T c corresponding to a critical coupling (3 C — 5.8. 
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m=0.100 ' — 
m-0.080 - a 
m=0.040 
m=0.020 
m=0.010 
m=0.000 




25 



20 



15 



10 



m=0.010 • 
m=0.020 
m=0.040 ♦ 
m=0.080 

m=0.100 




(b) 



6.2 6.4 P 6.6 

Figure 3: Fermion condensate (a) and fermionic susceptibility (b) calculated on 
lattices of size 8 3 x 4 for several values of the fermion mass. In Figure 3a we also 
show the result of an extrapolation to m = (see text). Also shown in Figure 3b 
are interpolating curves with the corresponding error band obtained with the help 
of the Ferrenberg-Swendsen reweighting method. 
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Figure 4: Chiral susceptibility versus 1 / y/m at three values of the gauge coupling 
in the intermediate phase, i.e. for < T < T c . 



mass 


0.02 


0.01 


/3c 


5.80 (3) 


5.77 (3) 



Table 4: Pseudo-critical couplings for the chiral transition determined from the 
location of peaks in the fermionic susceptibilities on lattices of size 8 3 x 4. 

The critical coupling, f3 c , may also be estimated from the location of the peaks in 
Xm, which become visible for small fermion masses. The corresponding pseudo- 
critical couplings extracted from a Ferrenberg-Swendsen interpolation are given in 
Table |4]. This is in accordance with our extrapolation of the fermion condensate. We 
thus obtain as an estimate for the critical couplings of the chiral and deconfmement 
transitions in the zero mass limit 

/3 d = 5.236 ±0.003 , 

(3 C = 5.79 ±0.05 . (15) 




The difference between these critical couplings may be used to estimate the ratio 
of the critical temperatures for deconfmement and chiral symmetry restoration in 
aQCD. To do so we use the 2-loop /^-function, Eq. [II]. This yields^ 



T c /T d « 7.7 ±2.1 . (16) 

f Using only the one-loop formula one obtains T c /Td = 11.9±2.8. Note that this is much smaller 
than the estimate given in Ref. j| for the case of SU(2), T c /T d = 175 ± 50. 
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We stress that so far no detailed investigation of the approach of lattice regularized 
aQCD to the continuum limit exists. If violations of asymptotic scaling turn out to 
be similar in magnitude to those found in the pure SU(3) gauge theory this ratio 
may well be underestimated by a factor of two. 



3.3 Latent heat and pressure close to 

In the previous subsections we have collected evidence for two separate phase transi- 
tions in aQCD. Let us now address in somewhat more detail the question how these 
transitions influence the equation of state of aQCD. The basic quantity we want 
to discuss here is the pressure as function of temperature. In numerical calculations 
one has direct access to the free energy density in units of T 4 , i.e. f/T 4 [T7|, [13. 
In the thermodynamic limit this is directly related to the pressure, p/T A = —f/T 4 . 
The free energy density is obtained from an integration over the difference of gluonic 
actions, {Sg), calculated at zero and non-zero temperature. The former only serves 
to normalize the pressure correctly, all significant changes in the pressure as function 
of temperature are, however, visible already in the expectation value of the gluonic 
action. In particular, the first order deconfinement phase transition also leads to a 
discontinuity (ASg), which directly is related to the latent heat at Td, 



T 4 T 4 da \N L 

Here {ASg) = {Sg)- — (Sg)+ denotes the difference of action expectation values 
calculated at T d in the confined (— ) and deconfined (+) phase, respectively. Sg 
is given by Eq. | and — ad(3/da is the /3-function of aQCD. The latter should be 
determined non-perturbatively at the critical coupling (3d- As a first estimate for the 
latent heat it is, however, sufficient to use the asymptotic 2-loop relation, Eq. [LI], 
between the gauge coupling (3 and the lattice spacing a. A similar approach has also 
been used to estimate the latent heat at the deconfinement transition in a SU{3) 
gauge theory |H| and we thus can compare our results with this calculation. 



In Figure | we show the plaquette expectation values which are related to the 
expectation value of the gluonic part of the action through 

<*> = 1 - <Sc> ' (18) 

Even on the rather small lattice used by us (8 3 x 4) the rapid change at Td is clearly 
visible. A comparison with results obtained in the pure gauge theory (m — > oo) on 
larger lattices shows, in fact, that the discontinuity in {Sg) is much larger in the 
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0.65 - 



0.60 



0.55 - 



0.50 




Figure 5: Normalized expectation values of the gluon action defined by Eq. 18 
(plaquette expectation value) for different fermion masses on lattices of size 8 3 x 4. 
Also shown is the corresponding result obtained in the pure gauge sector (m — > oo) 
on a 16 3 x 4 lattice. 



light fermion mass limit of aQCD. This also is evident from the probability distri- 
bution of plaquette values, which shows in all clear double peak structure. 
Again we have used the Ferrenberg-Swendsen reweighting to construct the plaquette 
distributions at (3 c {m). As an example we show in Figure |]the result for m = 0.04. 



0.05 



0.04 



0.03 



0.02 



0.01 



0.00 ' — 
0.48 



h(P) 



m=0.040 



0.50 



0.52 



0.54 



0.56 



0.58 



Figure 6: Distribution of the average plaquette values at j3 c (m = 0.04) on a 8 3 x 4 
lattice obtained by reweighting results at nearby /^-values. 



The difference in the location of the two peaks in this distribution gives an estimate 
of the discontinuity in the plaquette expectation value. The results for different 
fermion masses are summarized in Table ||]. A linear fit in the mass parameter yields 
the result at m — 0, which is given in the last column of Table |^. From this one 
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m 


0.10 


0.08 


0.04 


0.02 


0.00 


(AP) 


0.0220 (9) 


0.0227 (12) 


0.0335 (29) 


0.0308 (13) 


0.0336 (16) 



Table 5: Estimate of the discontinuity in plaquette expectation values at for 
several values of the fermion mass. The last column gives the result of a linear 
extrapolation to the zero fermion mass limit. 

obtains 

^ = 12.39 ±0.59 (19) 

as an estimate for the latent heat of the deconfmement transition. 

In order to compare with corresponding results in the pure gauge sector we should 
relate the latent to the energy density of a free gas of gluons and adjoint fermions. 
Taking into account the known cut-off dependence of a free Bose and Fermi gas on 
lattices with a small temporal extent N T = 4 |20| we find Ae/esB = 0.306 (15). This 
can directly be compared with the result in the pure gauge sector fll9"f , Ae/e SB = 
0.31 (3). We thus find that the latent heat per partonic degree of freedom is quite 
similar in the pure gauge theory and aQCD although in absolute units they differ 
quite a bit because 2-flavour aQCD has 2 • (iV 2 — 1)(1 + 2 • n/) = 80 degrees of 
freedom whereas the pure gauge theory only has 2 • (iV 2 — 1) = 16. 

Let us finally discuss the temperature dependence of the pressure, p/T 4 . It generally 
is obtained from the plaquette expectation values shown in Figure § by subtract- 
ing the zero temperature expectation values for each value of the gauge coupling 



j3 and then integrating this difference as a function of (3 [|T7| , jig] . Experience with 
similar calculations in the pure gauge sector shows that the zero temperature and 
finite temperature expectation values differ significantly only in the region where 
the finite temperature observable shows a rapid variation with j3. The zero tem- 
perature expectation values only provide the subtraction of a smooth background. 
Without going through the complete calculation of the pressure we thus can deduce 
its structure already from the finite temperature data shown in Figure [|. The ex- 
pectation values only show a rapid change at Tj. At the chiral transition point they 
are smooth. We thus find that also p/T 4 will increase rapidly only at and vary 
smoothly at T c . 

The chiral transition does not seem to leave a visible sign in bulk thermodynamic 
observables as it does not lead to a significant change in the number of light degrees 
of freedom at T<j as well as T c . The number of Goldstone modes, which should 
be present up to T c , is in the continuum limit about a factor eight smaller than 
the number of light adjoint fermions, which are liberated above However, our 
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calculation has been performed using the staggered fermion formulation, which only 
preserves a subgroup of the continuum SU(2rif) flavour symmetry. The actual 
number of light modes thus is reduced and the total number of nine light Goldstone 
particles will only be recovered in the continuum limit. This may reduce their 
influence on the bulk thermodynamics in our present analysis. Clearly this aspect 
of the thermodynamics deserves further studies closer to the continuum limit. 



4 Conclusions 



We have shown that a SU(3) gauge theory with fermions in the adjoint represen- 
tation leads to two distinct phase transitions at non-zero temperature. In the case 
of two Dirac fermions (simulated in the staggered fermion formulation) deconfine- 
ment occurs at a lower temperature than chiral symmetry restoration, T c /Td ~ 8. 
Unlike in QCD with fermions in the fundamental representation one thus has the 
possibility to examine the influence of both non-perturbative mechanisms on the 
thermodynamics separately. It will be interesting to analyze in the future in more 
detail the flavour dependence of T c /T d . In particular, it would be of interest to study 
both transitions in the supersymmetric limit, i.e. for rif = 1/2. 

We find that the deconfinement transition is first order for all values of the fermion 
mass and strongly influences the temperature dependence of bulk thermodynamic 
observables like the pressure as well as the screening of static charges in the funda- 
mental and in the adjoint representation. The latent heat per degree of freedom is 
similar to that found in a pure SU (3) gauge theory and also the screening of static 
fermion charges, which sets in immediately above T<j, becomes stronger than in the 
pure gauge sector and thus reflects the increasing number of partons (bosons and 
fermions) participating in the screening of external charges. 

The continuous chiral transition, on the other hand, does not seem to influence 
these observables significantly. We find, however, that chiral symmetry breaking 
in the deconfined region is significantly different from that in the confined region. 
Immediately above Td the chiral condensate starts showing a much stronger fermion 
mass dependence than below T d . The characteristic A/m-dependence, which leads to 
a diverging fermion susceptibility as expected from the analysis of three dimensional 
cr-models, is clearly visible for Td < T < T c . As we use in our present analysis 
the staggered fermion formulation it is likely that only a small subset of the nine 
Goldstone modes present in the continuum are light in our lattice calculation. A 
more detailed analysis of the Goldstone modes closer to the continuum limit would 
certainly be of interest in the future. 
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